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Abstract 

Errors are found in example problems from Henri Poincare's paper 
"Memoire sur les courbes definies par une equation differentielle." Exam- 
ples four and five from chapter seven and examples one, two, and three 
from chapter nine do not have the limit cycles at infinity predicted by 
Poincare. Instead they have fixed points at every point at infinity. In 
order to understand the errors made by Poincare, examples four and five 
are studied at length. Replacement equations for the fourth and fifth ex- 
amples are suggested based on the supposition that terms were omitted 
from Poincare's equations. 

1 Introduction 

In the late nineteenth and early twentieth century Henri Poincare began to study 
the quahtative aspects of systems of differential equations. This analysis was a 
breakthrough in the field because one no longer had to obtain a specific solution 
of the equation in question to understand its general behavior. This manner 
of analysis was introduced by Poincare in his paper "Memoire sur les courbes 
definies par une equation differentielle" |H1 1881]. Poincare worked primarily 
with systems of two variables and played a key role in identifying the existence 
of limit cycles with the Poincare-Bendixson theorem. Poincare also searched for 
a complete global analysis of a system of two variables; to do so he introduced 
analysis at infinity by means of the Poincare Sphere. These two aspects of his 
analysis join in a behavior known as a limit cycle at infinity. 

To illustrate the application of these techniques, Poincare presents example 
problems in chapters seven and nine. In chapter seven, Poincare presents five 
example problems. He states that examples three, four, and five have limit 
cycles at infinity. (Examples three, four and five throughout the paper will be 
referred to as P3, P4, and P5 respectively.) However, when one checks Poincare's 
assertion, one finds that examples P4 and P5, in fact, do not have limit cycles 
at infinity. The errors remains in the collected works [7| and the error in P4 
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remains in a modern work that includes translations of examples P3 and P4. 
p. 155]. 

In chapter nine, Poincare presents three additional example problems. He 
states that each of these examples have limit cycles at infinity. However, he 
makes the same error analyzing these equations as he did for P4 and P5. 

Throughout this paper, in an effort to understand why Poincare made these 
errors, we focus on studying examples P4 and P5 from chapter seven. In the 
conclusion we will attempt to apply what we have learned about the errors in 
P4 and P5 to the examples from chapter nine. 

Because the analysis at infinity is algebraically simple for each of these prob- 
lems, one becomes perplexed at this finding. Surely Poincare could not have 
made such an error. With this belief, one might suppose that there was a simple 
omission of some sort. If one were to find equations that have the same behavior 
as Poincare predicts for P4 and P5, and that, were one to drop certain terms, 
could become his printed equations, this would serve as a plausible explana- 
tion. We look at P4 in the context of P3, and consider a plausible modification. 
Such an analysis is encouraged by the strong geometric similarities between the 
two systems. The modification restores the limit cycle at infinity, and hence 
supports the hypothesis of an error of omission, becoming a key candidate for 
the intended fourth example. We will call this example R4. Example P5 does 
not fit in the same genre of geometric system and hence no such modification 
suggests itself. 

When one investigates the uniqueness of these "solutions" to the problem of 
finding modifications of Poincare's equation that match his analysis for P4, one 
can find a classification of all modifications of the original equation, of a certain 
type, that result in the correct behavior. There exist a great number of such 
modifications; however, this remarkable non-uniqueness in "solutions" does not 
cause skepticism that R4 is the intended equation for the fourth example. It is 
seen that the modification to obtain R4 is the simplest of the modifications found 
under the classification, and hence, through application of the philosophical 
principle known as Occam's razor, R4 remains the most plausible equation for 
the fourth example. 

Based on the success of the technique applied to P4, a similar technique is 
applied to P5 to try to determine a possible modification of P5 which yields 
the correct behavior. This technique provides a number of modifications which 
agree with Poincare's analysis for P5 thus illustrating the power of the technique. 
One of these, which we shall call R5, is chosen as a suggested replacement to 
P5. (No modifications are suggested for the examples from chapter nine.) 

To aid the reader with the analysis of differential equations at infinity, a sec- 
tion is included presenting the standard techniques for analysis of the behavior 
at infinity for a polynomial system in rectangular coordinates. A more geomet- 
ric analysis is also presented that is quite valuable in many cases in which such 
a system simplifies in polar coordinates. 
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2 Differential Equations at Infinity 



Differential equations that are defined on the plane can be analyzed at infinity 
by extending the plane to a representation of the plane that includes the "points 
at infinity." This analysis began with the work of Poincare and Bendixson |H|[2]- 
Throughout this paper we use Poincare's method of analysis. More modern and 
readable accounts on the Poincare's method are available in books by Perko [S], 
Lefschetz [3] , and Minorsky (5j . This section begins with a brief summary of the 
standard results as in Perko [BJ pp. 264-268]. A theorem is presented, consistent 
with the standard techniques, based on polar coordinates, that demonstrates 
the geometric nature of these systems at infinity. This theorem often makes it 
possible for one to determine by inspection whether a system has a limit cycle 
at infinity. 

Through the use of projective geometry, the complete behavior of a two 
dimensional system of differential equations can be seen in its behavior on a 
sphere of finite radius known as the Poincare sphere. To do this, one places 
the phase plane tangent to the sphere and makes correspond points on the 
plane with the points on the sphere by central projection (a point on the plane 
corresponds with an antipodal pair on the sphere.) One must note that the line 
intersects the sphere at two points; to remove this non-uniqueness, antipodal 
points are identified on the Poincare sphere. The points that were at infinity on 
the original plane become points on the equator of the sphere. 

From the behavior on the Poincare sphere, one can construct the "global 
phase portrait." To consider the global phase portrait of a system, one projects 
the trajectories from the upper hemisphere orthogonally down onto the plane 
that goes through the center of the sphere and is parallel to the original plane. 
In this way the complete behavior on the original plane becomes the behavior 
on the finite disk of this new plane. Points at infinity become the points at 
which the sphere intersects the plane - the boundary of the disk. 

The first analysis that one would do for a system in the finite plane is to find 
the location of the fixed points. With this as motivation, the primary problem 
of analyzing the behavior of a system at infinity is to determine the location of 
fixed points, if there are any. One considers the system: 



where P and Q are polynomials of degree m in x and y. Denote by P„i and Qm 
the homogeneous polynomials consisting of terms of degree exactly to. 

Theorem 1 The critical points at infinity for the system above occur at the 
points (X,Y,0) on the equator of the Poincare Sphere where + = 1 and 



dx 

dy 
dt 



XQ„,{XX)-YP,n{X,Y) = 
or equivalently at the polar angles 6 and 6 + t: which are solutions of 
Gm+i{Q) = cos6 Q,n{cos9, sin9) — sin6 P„i{cos6 , sin6) = 0. 



(2.1) 



(2.2) 
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This equation has at most m + I pairs of roots 6 and 9 + n unless Gm+i{0) 
is identically zero. If Gm+i{0) is not identically zero, then the flow on the 
equator of the Poincare sphere is counter-clockwise at points corresponding to 
polar angles 6 where Gm+i{d) > and it is clockwise at points corresponding to 
polar angles where Gm+i{0) < 0. 

See Perko 0, section 3.10, Theorem 1. 

An important consequence of this Theorem 1 is that any polynomial system 
in rectangular coordinates can be extended to the Poincare sphere. 

One should be aware that in the proof of Theorem 1, a new time scale is 
defined to study the trajectories in the proximity of the equator. There are two 
opposing conventions used in the degenerate case where (|2.1() is identically zero. 
The most common one, used by Perko jBj and Minorsky (Sj is to re-parameterize 
time in such a way that the equator is forced to consist of trajectories and fixed 
points. The second convention, followed by Lefschetz ^ and Poincare jH| uses 
a different time scale, allowing trajectories to cross the equator (implying that 
trajectories in one of the hemispheres will flow in a direction inconsistent with 
the flow on the plane.) Because the former method is more common, we will 
adhere to it, even though Poincare used the latter. (A brief comment will be 
made in the analysis of Example P4 about why this does not invalidate our 
results.) 

If there are no fixed points at infinity there is a cycle at infinity. If trajectories 
in the proximity of this cycle at infinity approach (or recede from) it, it is a limit 
cycle at infinity. A corollary immediately follows: 

Corollary 1 Let r denote radial distance in polar coordinates. If the system 
considered above has no fixed points at infinity and dr/dt ^ for r > R, for 
some R, then there exists a limit cycle at infinity. 

While this theorem makes it relatively easy to determine the behavior of a 
system at infinity, there is a more geometric analysis that for certain systems 
will make the behavior at infinity particularly easy to determine. 

Theorem 2 Consider a polynomial system 

^ = P(:r,y), f = (2-3) 

with expressions for dr/dt and dd / dt of order I and J in r as r ^ oo, respec- 
tively. Let k — I— J. Then, if k>2, the equator of the Poincare sphere consists 
entirely of fixed points. If k < 1 then: if Gm+i{()) 7^ for all 6 there is a cycle 
at infinity. Furthermore, if dr/dt satisfies the conditions of the above corollary, 
it is a limit cycle. Otherwise, the equator of the Poincare sphere has finitely 
many fixed points located at 9 such that Gm+i{9) = 0, as in Theorem L 

Proof of Theorem 2: 

Suppose that P and Q are polynomials of degree m in a; and y. We can express 
this system in polar coordinates, for r 0, as follows, 

dr 

— = (cos9Pq + sin9Q{)) + r(cos9Pi(cos9, sin9) + sin9Qi(cos9, sin9))-\- 
dt 
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• • • + {cos9 Pm{cos9 , sinO) + sin9Qm{cos6, sinO)). 

do 

— — r^^ (cosOQo — sinOPo) + {cos9Qi(cos6, sinO) — sin9Pi(cos6, sin6))+ 
dt 

• • • + r"^^^ {cos9Q„^{cos9 , sin9) — sin9Pm{cos9, sin9)). 

Throughout the remainder of the proof we require r ^ 0. The foUowing defini- 
tions wiU greatly simphfy notation: 

770(6*) = {cos9Po+sin9Qo), 

rii{9) — {cos9 Pj {cos 9, sin9) + sin9Q I {cos9, sin9)), 
^_i(6») = {cos9Qo- sin9Po), 



and 



Cj(^) = {cos9Qj+i(cosO,sin9) — sin9Pj+i{cos9,sin9)), 

where I and J are the highest degree terms, in r, in the equations for dr/dt and 
dO/dt respectively. With this notation we obtain, 

J - ^,{9) + --- + r'm{0), 

^ = r-'^-m+U0) + --- + r'U9). 

By the definition of k we have k = I — J. If wc put this system in differential 
form we obtain 

(r-i^-i W + • • • + r'O{0))dr - + • • • + rS/(^))d^ = 0. (2.4) 

To determine the behavior at infinity, one projects a differential equation 
onto the Poincare sphere. To understand the details of the behavior at infinity, 
project the upper hemisphere of the Poincare sphere onto the cylinder of radius 
1 with axis orthogonal to the phase plane that is tangent to the Poincare sphere 
at the equator. (We can restrict our attention to the upper hemisphere because 
we follow the standard convention, which does not allow trajectories to cross the 
equator.) The "equator" of the cylinder is the part of the cylinder that touches 
the sphere. By geometric analysis of the projection from the plane directly to 
the cylinder, s = 1/r and correspondingly dr = ^ds. This projection leaves 9 
unchanged. 

Projecting (|2.4II from the plane directly onto the cylinder gives: 



(s^^m + ■■■ + ^ij{9))i^ds) - {m{9) + ■■■ + ^Vii0))d9 = 0, 
which simplifies to 

{s^-i{9) + ■■■ + ^Omds + s\f^o{0) + ■■■ + ^Vimd9 = 0. (2.5) 
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Case 1: k>2 

In this case we have / — J > 2, so / — 2 > J. Wc multiply (|2.5() by s^~'^ to 
clear the denominators obtaining: 

{s'-^^-i{0) + ■■■ + s'-'-^tj{0))ds + {s'rjoiO) + ■■■+ m{0))de = (2.6) 

However, this equation would indicate trajectories crossing the equator (from 
one side or the other) for all points [O^Q] where rji{9) ^ 0. Because we wish 
to maintain the same convention as in Theorem 1, the equator is required to 
consist only of fixed points and trajectories moving along the equator. As a 
result, the trajectories crossing the equator would lead to a violation of the 
uniqueness of solutions to (|2.6|l at these points (0,0). To resolve this, we must 
multiply H2.6|l by an additional value of s, a reparameterization of time that 
causes these trajectories to slow down so that they do not cross the equator. 
We obtain: 

(/e-i(^) + • • • + s'-^-'^j{e))ds + {s'+^Vo{0) + ■■■ + sT^i{9))d9 = (2.7) 
With regards to this new time scale, r, (|2.7|) can be expressed as the system: 

^ = (/^_i(0) + ••• + , s^-^-iO(e)). (2.8) 
Because fc = / — J > 2, at s = 0, we have: 

dr dr 

We conclude that every point on the equator is a fixed point. 
Case 2: fc < 1 

We have / - J < 1 and J < / - 1. We multiply (|2.5(l by s'-' to clear the 
denominators. In this case, no trajectories approach the equator in a finite time, 
hence there is not need for the re-parameterization done in case 1. 

We obtain: 

(s^+ie-i(^) + - • ■ + s^j-i{9)+O{0))ds+{s-'+^M0) + - ■ ■ + s'-^'-'^i^i{9))d9 - 0. 

Expressing this in the form of a system, with respect to the new time scale r, 
we obtain: 

di 

^ = -s^-''-'\s%,{9) + ---+r^j{9)l 

^ = is'+^^_,[9) + --- + ue)). 

So, on the equator we have s — Q and 

dr dr 
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So, the fixed points at infinity are for 9 such that ^j{0) = 0. It is easy to see 
that, in this case, ^j{9) = Gm+i{9), as defined in Theorem 1. Hence, the fixed 
points at infinity are at 9 such that 

Gm+i(^) ~ cos9Qrn,{cos9, sin9) — sin9Pm{cos9, sin9) — 0. 

We conclude that if Gm+i{9) ^ for all 9 then there is a cycle at infinity, and 
that, if nearby trajectories approach that cycle (or recede from it), as in the 
corollary to Theorem 1, then it is a limit cycle. 
This concludes a proof of Theorem 2. 

One should notice that in the proof of this theorem one finds that for the case 
k < 1 one has Gm+ii9) = ^j{9), the highest order term in r of the d9/dt equa- 
tion. This makes this theorem particularly useful - one can often tell whether 
a system has a limit cycle at infinity by merely expressing the system in polar 
coordinates. 

One should also note that in the proof of case 1, if fc > 2, equation (|2.8|l 
gives that a trajectory approaching (or receding from) any point on the equator 
does so orthogonally to the equator. 



3 The Poincare Examples 

To enable the reader to understand Poincare's fourth and fifth examples, three 
of Poincare's examples from chapter seven of his paper jj) pp 274-281] are 
presented. Poincare's third example, P3, works as an introduction to his fourth 
and fifth examples and serves as an example of a system that has a limit cycle 
at infinity, the behavior that Poincare claims for P4 and P5. Examples P4 and 
P5 are then presented with demonstration of the error in his analysis at infinity. 
At the end of this section, the error in Poincare's analysis of his examples from 
chapter nine is briefly discussed. 

Example P3 

Poincare considers the equation: 

dx dy 



x{x'^ + y'^ — I) — y^x"^ + + 1) yix^ + y"^ — 1) + x{x'^ + j/^ + 1) 

This differential form of the equation is equivalent to the more familiar form of 
a system as noted in (Bj p. 266]. 

^ = xix^ +y^ -l)-yix^ +y^ + 1), 

^ = y{x^ +y^ -l) + xix^ +y^ + 1). 

The geometric aspects of this system become far more obvious in polar form: 
dr / 2 IN 



'dt 



+ 1. 
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Poincare presents the key qualitative features of this system: 

• The origin: (x, y) = (0, 0) is a fixed point. One has a stable spiral at the 
origin. 

• A repelling limit cycle at radius r = 1 centered about the origin. 

• An attracting limit cycle at infinity. 

All of the features in the finite plane are easily established, so only the features 
of this system at infinity are analyzed in detail. Poincare establishes that there 
is a limit cycle at infinity because he claims that there are no fixed points at 
infinity: "II n'y a aucun point singulier sur I'equateur, qui est une caracteristique 
et qui est par consequent un cycle limite" '5', p. 278]. He omits this calculation, 
so we verify it here. 

We apply Theorem 2 to show that there are no fixed points at infinity. From 
the expression for the P3 in polar coordinates, we find / = 3 and J — 2. Further, 
^2(^') = 1 for all 9, demonstrating the existence of a cycle at infinity. Because 
dr/dt > when r > 2, this cycle is a limit cycle. 

Example P4 

Poincare's fourth example can be considered a more complicated version of 
his third example. Poincare writes the equation as: 

dx 

x{x'^ +y'^ - l)(a;2 + y2 _ g) _ y(a;2 + y2 _ 2a; - 8) ~ 

dy 

y{x'^ +y'^ ~ + y2 _ 9) + x(^x2 + y2 _ 2a; - 8) ' 

which he transforms into the system: 

^ = x{x^ -l){x^ +y^ -9)~y{x^ +y^ -2x~8), 

§ = yix^+y^-l){x^+y^~9)+xix^+y^-2x-8). 
dt 

Poincare states the key features of the system: 

• The origin: (a;, y) = (0, 0) is a fixed point. One has an unstable spiral at 
the origin. 

• A: {x,y) = (1/2, V35/2), and B: {x,y) = (1/2, -\/35/2) are fixed points: 
an unstable node and a saddle respectively, at the intersection of the circles 
defined by + - 9 = and a;^ + - 2a; - 8 = 0. 

• An attracting limit cycle of radius r = 1 centered about the origin. 

• Two heteroclinic orbits connecting fixed points A to fixed point B along 
the circle + _ g = q. 

• An attracting limit cycle at infinity. 
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Figure 3.1: Sketch of the global phase portrait corresponding to Poincare's 
analysis of example P4. Circles Ci,C2, and C3, which are referred to later in 
this paper, are labeled. (The dashed line for C3 indicates that it is a nuUcline 
and not a trajectory.) 

All of the stated features in the finite plane can be easily established. (A 
similar analysis is applied to a slightly more complicated system in the following 
section number five.) The geometric nature of the system becomes clear in polar 
coordinates. The system becomes: 



The phase portrait corresponding to Poincare's analysis is included in Figure 3.1. 
He claims that, as in P3, there are no fixed points at infinity, and hence a limit 
cycle. However, when one does the calculations to check Poincare's assertion, 
one finds that this is not the case. 

To check the assertion, one applies Theorem 2 to the system. By looking 
at the the expression for P4 in polar coordinates, one finds that 1 = 5 and 
J = 2, hence k = 3 and every point at infinity is a fixed point. There is no limit 



dr 

'di 

dO 

'dt 



r{r^ - l)(r2 - 9), 



- 2rcos9- 8. 
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Figure 3.2: The correct global phase portrait corresponding to Example P4. 
(Figures 3.2 and 4.1 were calculated numerically using a Runge Kutta method.) 

cycle at infinity because every point on the equator is a fixed point. (If one had 
used the second convention, which was used by Poincare, one would have found 
trajectories crossing the equator. With this convention, the trajectories crossing 
the equator eliminate the possibility of a cycle on the equator. All reference to 
this second convention will be dropped for all of the following examples.) The 
correct global phase portrait for the system is presented in Figure 3.2. 
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Example P5 

Poincare's fifth example is not as clear of an extension of his third and 
fourth examples as it does not quite fit the same "genre" of system. The main 
difference is that whereas P3 and P4 arc more naturally understood in a polar 
coordinate system, P5 is more easily understood in a coordinate system based 
on lemniscates. Example P5 is also more complicated because it depends upon 
the parameter c. He considers the equation: 

dx 

a;(2x2 + 2j/2 + l)((a:2 + y2)2 + ^2 _ y2 _ c) - y(2x'^ + 2y2 _ 1) = 

dy 

y{2x^ + 2y2 - + y2)2 + 3,2 _ ^2 _ ^) _|_ a;(2a;2 + 2y^ + 1)' 

Poincare asserts that P5 has the following behavior, which depends on the 
parameter c. 

1. For c < —1/4, the system has: 

• A saddle point at the origin. 

• Two unstable spirals at (0, ±l/\/2). 

• An attracting limit cycle at infinity. 

2. For — 1/4 < c < 0, the system has: 

• A saddle point at the origin. 

• Two stable spirals at (0, ±l/\/2). 

• Two repelling limit cycles surrounding each of the two above spirals 
respectively. 

• An attracting limit cycle at infinity. 

3. For c = 0, the system has: 

• A saddle point at the origin. 

• Two stable spirals at (0,±l/\/2). 

• Two homoclinic orbits which branch from the saddle and surround 
each of the spirals respectively. 

• An attracting limit cycle at infinity. 

4. For c > 0, the system has: 

• A saddle point at the origin. 

• Two stable spirals at (0,±l/\/2). 

• One single repelling limit cycle centered at the origin which surrounds 

all three fixed points. 

• An attracting limit cycle at infinity. 
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Case 1 




Figure 3.3: Global phase portraits corresponding to Poincare's analysis for the 
four different cases of P5. 
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The four qualitatively different phase portraits are sketched in figure 3.3. As 
in the example P4, all of these features in the finite plane can be easily verified. 
Because Poincare does not include his analysis of the behavior at infinity we 
check it here. Example P5 does not simplify in polar coordinates, so it is easiest 
to use Theorem 1 to determine the behavior at infinity. The highest order terms 
are of degree seven so we have: 

Qr = 2Y{X^+Y^)^. 

So to find the fixed points at infinity, according to Theorem 1, one must solve 
the following system: 

XQr ~ YP-j = 2XY{X'^ + Y'^f - 2YX{X'^ + Y'^f ^ 0, 

X^ + Y^ = 1. 

The first equation is satisfied for all X,Y, so every point at infinity is a fixed 
point. Hence there can be no limit cycle at infinity. (This error was first 
discovered by Professor Lawrence Perko of Northern Arizona University when 
he read a preliminary draft of this paper.) 

Examples from chapter nine 

Poincare presents three additional examples in chapter nine jSl p. 292-296] of 
his paper to illustrate the more general techniques that he developed in chapter 
eight. Poincare claims that each of these examples has a limit cycle at infinity. 
However, one can easily check his claim using Theorem 1 and Theorem 2, as 
was done above for P3, P4, and P5, to see that each of these examples has 
fixed points at every point on the equator. Because it is easy to check this, the 
verification is left to the reader. Furthermore, throughout this paper, we focus 
on a discussion of examples P4 and P5. 

The fact that these errors were so easily detected immediately makes one 
wonder whether Poincare used a similar form of analysis for finding fixed points 
at infinity as was presented in Section 2. When one looks at his first and second 
examples P, p. 274-278] one finds that he uses the equations from Theorem 1 to 
find fixed points at infinity. Because the calculations for finding fixed points at 
infinity for each of these examples are so simple (using Theorem 1 or Theorem 
2) one is led to speculate as to the cause of Poincare's error. 

4 Error of Omission 

The facts that Poincare's errors in finding the fixed points at infinity for his 
examples P4 and P5 were so easily detected, and that he used the same type 
of criterion that is used here, suggest that, perhaps, his analysis was of differ- 
ent equations than the ones printed in the paper. Were there to be an error 
of omission, terms which would lead to an unsolvable system for finding the 
fixed points at infinity could be missing which would not change the qualitative 
nature of his system on the finite plane, where his analysis matched the printed 
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equation. If this were the case, there most probably was an error of omission in 
the production of his manuscript, not in the mathematics. After lengthy inves- 
tigation, a system which could very plausibly become the system printed in the 
paper as P4. were one to omit certain terms, becomes convincing candidates for 
the intended fourth example. 

In terms of the qualitative features predicted, Poincare's fourth example ap- 
pears to be an augmented version of his third example; however, the equation 
listed is not an algebraically augmented version of the equation in his third 
example. Algebraically, when in polar coordinates, the dr/dt equation is aug- 
mented with another nullcline at r = 3. To create the new fixed points A and 
B, the dO/dt equation has a new nullcline which crosses the dr/dt nullcline at 
r = 3. However, the d6/dt equation is missing the factor (r^ -|- 1), which would 
not create any nuUclines, but which was necessary for P3 to have a limit cycle 
at infinity (see calculation.) A plausible attempt at finding Poincare's actual 
equation is to include this term in the polar version of example P4. This at- 
tempt would fit with the intuition that P4 is merely an augmented version of 
P3. (Example P5 is not considered "in the context of examples P3 and P4" 
because the limit sets for P5 are based on lemniscates instead of circles.) 

Example R4 

The system becomes: 

I = r(r^-l)(r^-9), 

^ = (r2-2rcos(0)-8)(r2 + l). 
dt 

In Cartesian coordinates the system becomes: 

^ = xix" -l){x' -9)-y{x' +y' -2x-?>){x^ + + 1), 

^ = y(a;2 + y2-l)(a;2+t/2-9) + a;(a;2+y2-2a;-8)(a;2+t/2 + l). 

This new version of Poincare's fourth example has all of the qualitative features 

that his printed system has on the finite plane, and has the additional feature 
of a limit cycle, instead of fixed points, at infinity. Hence this system exactly 
matches Poincare's analysis of P4. The global phase portrait corresponding to 
example R4 is included in figure 4.1. Verification that this system has the same 
qualitative features as Poincare claims for example P4 will not be presented here 
because it is a result of Theorem 3 in the following section. Because this system 
algebraically follows from Poincare's third example, according to the hypothesis 
that Poincare merely augmented his third example to obtain his fourth example, 
this compels the author to believe that this was Poincare's intended system. 

5 A Method of Classification 

The results of the last section raise the question: how can one change the 
algebraic statement of a differential equation and have the qualitative aspects 
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Figure 4.1: The global phase portrait corresponding to example R4. 
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in the finite plane remain the same, while changing the behavior at infinity? 
Rather than using ad hoc methods for changing an equation, can there be a 
more systematic way of finding a certain modification of an equation with the 
desired conseqiicnces? 

To begin to answer these questions this section contains a systematic way 
for finding which, among a certain class of changes, cause the equation listed in 
Poincarc's paper as P4 to be qualitatively equivalent to his analysis listed for 
P4. Based on the success of this technique with P4, it is applied to P5 leading 
to a recommended modification of P5 which we will call R5. These results not 
only serve as a demonstration of such techniques, but also as a statement about 
the remarkable non-uniqueness of a equations with this specific phase portrait. 

This is the modification of P4: 

Theorem 3 (Classification Theorem for P4) Let W{x, y) be a polynomial 
in X and y of degree N, and let Wjv(a;,y) be the homogeneous portion of W 

which is of degree N . If 

1. W{x,y) > 0, for all {x,y) e R^, 

2. N>2, and 

3. WN{x,y) is positive definite. 

then: the phase portrait of the system given by: 

x{x^ +y^ - l){x^ + y2 - 9) - y{x^ + y^ - 2x - 8)W{x, y), 

yix"^ + 2/2 _ i)(a;2 + y-^-Q)+ x{x^ +y^-2x- 8)W{x, y). 

gives a phase portrait qualitatively equivalent to the analysis Poincare gives in 
example P4- 

First, one must verify that the system has the same number, location, and 
type of fixed points. To do this, first convert the system to polar coordinates, 
obtaining: 

r{r^ - l)(r2 -9), 
(r^ — 2rcos6 — 8)W{rcos6, rsinO). 

Because the new factor in the d6/dt equation is never zero, it does not create 

any new nuUclines. Hence, because there are no new nuUclines, only the fixed 
points that existed in the equation listed in the paper {A,B, and O) exist in 
these modified examples. 

To verify that these fixed points have the same local behavior as presented 
by Poincare one must calculate the trace, r, and determinant. A, of the Jaco- 
bian matrix of the system (in rectangular coordinates) symbolically at arbitrary 



dx 

Itt 
dy 

dt 



dr 
dt 
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{xo,yo)- To simplify the calculation consider the following definitions: 

Ci{x,y) = + 
C2{x,y) = + 
Cz{x,y) = x"^ + y"^ - 2x - 8. 

With these definitions our system becomes: 
dx 

— = xCi{x,y)C2{x,y) -yC3{x,y)W{x,y), 
dy 

— = yCi{x,y)C2{x,y) + xC3{x,y)W{x,y). 

The expressions are calculated as follows (the arguments of circles Ci, C2, and 
C3 are dropped): 

_ dP_ dQ 
dx dy ' 
^ _ dPdQ _dPdQ 
dx dy dy dx ' 

One obtains: 

T = 2CiC2 + 2{x^+y^){Ci+C2) + 2yW{x,y) + 

dW{x,y) dW{x,y) ^ 
^^^"^ dy y dx 

(C1C2 + 2x^C, + C2) - y{2x - 2)W{x, y) - yC^^^^)- 

(C1C2 + 2y\C^ + C2) + 2xyW{x, y) + xC^^^^)- 
{2xy{C^+C2) - C3W{x,y) - 2y^W{x,y) - yC^^^^)- 
{2xy{Ci + C2) + C^Wix, y) + x{2x - 2)W{x, y) + xC^^}^^). 

When one evaluates these quantities at one of the fixed points, to check 
the Unearization, one is confronted with the unknown function W{x,y). One 
must make restrictions on W{x, y) so that the linearizations satisfy Poincare's 
analysis for P4. 

At fixed point O all terms multiplied by a; or y vanish, so one obtains: r = 18, 
A = 81 + 64W(0,0)^ and t2-4A = -256W(0,0)^ Because r > and because 

— 4A < (because W{Q, 0) > 0) we have an unstable spiral. 

At fixed point B we have A = -144^351^(1/2, -a/35/2). Once again, 
because W{x, y) > for all {x, y), we have A < 0, hence a saddle. 

At fixed point A we have r = 144 + 7351^(1/2, \/35/2) and 
A = 144V35VF(l/2, \/35/2). Because W{x, y) > wc have t > 0. Furthermore, 

- 4A = (144 - 7351^(1/2, 735/2))^ > 0. 

Because — 4A > and r > 0, we find that fixed point A is an unstable node. 

Hence for any W satisfying the conditions that W{x,y) > for all {x,y) 
then the fixed points of the modified system have the same behavior as listed 
in Poincare's analysis of P4. 



A = 
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It is easy to verify that there is a stable Hmit cycle at r = 1. One must 
notice that dr/dt = r{r'^ - l)(r^ - 9) = for r = 1, dr/dt > for < r < 1, 
dr/dt < for 1< r < 3, and that dO/dt = (r^ - 2rcose-8)W{rcose, rsinO) < 
for < r < 2. Hence r = 1 is a stable limit cycle. 

To establish the existence of the heteroclinic orbits, notice that for r = 3, 
dr/dt = 0. On this circle of radius three, to the left of A and B, dO/dt > 0, 
and to the right of A and B, dO/dt < 0. The sign of dd/dt remains the same 
in these regions as in the printed equation because W{x,y) > for all {x,y) in 
R^. Hence there are trajectories which maintain r = 3, connecting ^ to S both 
directions along the circle. 

Finally we apply Theorem 2 to establish the existence of a limit cycle at 
infinity. We can express W{x, y) in polar coordinates in the following way: 

W{rcos6, rsinO) = Wo{cos6, sinO) + rW\{cos6, sinO) H h r^WN{cos6, sinO). 

With Wn{cos6, sinO) > for all 9 because Wn{x, y) is positive definite. Hence, 
using the notation of Theorem 2, / = 5 and J ^ 2 + N, where the degree of 
W{x,y), N > 2. Furthermore, ^j{9) = Wn{cos0, sin9) ^ for ah 9. Hence, 
there is a cycle at infinity. This cycle is a limit cycle because for r > A, dr/dt > 0. 
This concludes the proof of Theorem 3. 

Based on the success of the above method one may wish to see how well it 
works for another system-we use P5. 

Theorem 4 (Classification Theorem for P5) Let W{x, y) he a polynomial 
in X and y of degree N, and let WN{x,y) be the homogeneous portion of W 
which is of degree N . If 

1. W{x,y) > 0, for all {x,y) e R^, 

2. N>4, and 

3. WN{x,y) is positive definite, 

then: the phase portrait of the system, given by: 
dx 

— = x{2x^ + 2y^ + l){{x^+yY+x^^y^-c)-y{2x^ + 2y^-l)W{x,y), 

5 = y{2x^ + 2y^-l){{x^+yY + x'-y^-c) + x{2x^ + 2y^ + l)W{x,y), 
dt 

gives a phase portrait qualitatively equivalent to the analysis Poincare gives in 

example P5. 

We follow the same technique as in the previous theorem; first we verify 
that the modification results in the same fixed point behavior. Clearly, W{x, y) 
is never zero, so it docs not lead to any new fixed points. Hence we must 
only check that the existing fixed points have the same behavior. To make the 
analysis easier, we make the following definitions: 

A{x,y) = x{2x'' + 2y^ + l), 
B{x,y) = 2/(2x2 + 22/2-1), 
C{x,y) = {x'+yy + x^-y^-c. 
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With these definitions the system becomes: 
— = A{x,y)C{x,y)-B{x,y)W{x,y), 

^ = Bix,y)C{x,y)+Aix,y)Wix,y). 

At the fixed points A{x, y) = B{x, y) = and x = so in the expressions for 
the trace and the determinant we have: 

T = 8y^C{x,y), 

A = i2y' + l){6y'-l){{C{x,y)f + {W{x,y)f) 

At (0,0) we have A = — (c^ + W{0,0)^) which is negative for all {x,y) 
because W{0, 0) > 0. Hence wc have a saddle at the origin, independent of the 
parameter value c and the modification W . 

At (0,±l/\/2) wc have t = 4C(0,±1/V2) = -4(1/4 + c) and 
A = 4C(0, ±l/V2)^+4W(0, ±l/\/2)2. Hence, -4A = -16W{0, ±1/V2f < 
resulting in spirals for any choice of W. The stability depends upon t; here 
T > for c < —1/4 and r < for c > — 1/4. This matches Poincare's assertion 
of imstablc spirals for c < —1/4 and stable spirals for c > —1/4. Hence, as long 
as W{x, y) > for all {x,y) the behavior of the fixed points is unchanged by 
W. 

Next, we demonstrate that C(x, y) = is a conserved quantity. This results 
in the limit cycle(s) of cases 2 and 4 and the homoclinic orbits of case 3. Consider 
the rate of change in C{x,y) with respect to time: 

dCjx, y) _ dCjx, y) dx dC{x,y) dy 
dt dx dt dy dt 

= ^-^^{AC-BW) + ^-^^{BC + AW). 
dx oy 

Now one makes the observation that ^ = 2A and ^ = 2B. So that: 

^ = 2{A^+B^)C. 

So, clearly C{x,y) = is a conserved quantity. At this point we sec that the 
homoclinic orbits for case 4 have been verified because for c = the algebraic 
curve C{x, y) = is a leminescate based at the origin. Further, for C < 0, we 
have ^ < and for C > we have ^ > 0. Based on these results, one may 
find the necessary "trapping regions" to prove the existence of the limit cycles 
using the Poincare-Bendixson Theorem. 

Finally, Poincarc asserted that for all of these values of c there is a a limit 
cycle at infinity. To verify this, we determine whether there are fixed points at 
infinity by using Theorem 1. Using that N > 4:, we obtain: 

XQ^-YP„, = 2{X^ + Y^fWN{X,Y) = 0, 
X'^ + Y^ = 1. 
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which is clearly inconsistent because Wn{X,Y) is positive definite, by hypoth- 
esis. 

Now that it has been verified that there are no fixed points at infinity, to 
prove that there is a Hmit cycle at infinity it must be shown that trajectories near 
infinity do, in fact, approach the cycle at infinity. In this case we must apply 
slightly different techniques than Corollary 1 because dr/dt is not single signed 
for all values of 0. However, wc only must show that the system's behavior near 
infinity does approach the cycle at infinity. To do this we note that dC/dt = 
2{J^+B'^)C > for all points in the plane which are outside of the limit cycle(s) 
or homoclinic orbits defined by the algebraic equation C{x, y) = 0. Based on 
this, one can rule out other limiting sets in the proximity of infinity. Because 
there are no other limiting sets within a certain distance from the equator, the 
trajectories near the equator must approach the cycle at the equator-the limit 
cycle at infinity. 

This concludes proof of Theorem 4. 

These theorems demonstrate the incredible algebraic non-uniqueness of the 
equations which have the qualitative features outlined by Poincare for P4 and 
P5. This non-uniqueness might be of concern in speculating which modifi- 
cations are the right ones. However, were one to speculate, without insight 
into Poincare's other examples, as to which of these possible modifications 
to P4 might be his equation, one might first try the most simple of all W, 
W{x,y) = (x^ + + 1). This gives example R4 exactly. It is interesting that 
through choice by simplicity (application of Occam's Razor) one obtains the 
same result which was found through analysis of the context of the problem. 
This interesting philosophical observation continues to persuade the author that 
R4 really was the system intended by Poincare for P4. 

Based on all of the above success, the author would like to suggest an ade- 
quate system to replace P5: 

Example R5 

The system becomes: 

J = x{2x'' + 2y' + l){{x' + y^f +x^-y^-c)- y{2x^ + 2y^ - + / + 1), 

5 = y{2x^ + 2y^ - \){{x^ + y^f +x^-y^-c)+ x{2x^ + 2y^ + l){x^ + y^ + l). 
dt 

The fact that R5 matches Poincare's analysis of P5 is a clear application of 
Theorem 4. 

6 Discussion 

The reader may wonder why these problems were investigated. This section 
briefiy answers this question and attempts to bring closure to the problem. 

Verification of Poincare's assertions about his fourth example were given 
as part of the take home project for a differential equations class. We were 
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given the English translation as a guide . One aspect of the assignment was 
to attempt, through the use of a computer, to obtain an actual global phase 
portrait. When this was done, the trajectories were found to not cycle around 
near the equator, but to come in virtually orthogonal to it! (See figure 3.2.) It 
is by this means that the error was first suspected. Much speculation took place 
before all of the analysis of P4 was done. 

After a complete analysis of P4 came about, the error in example P5 was 
discovered by L. Perko. With example P5 in need of a modification and with 
faith in the classification technique used for P4 (despite the fact that it only 
considers one type of modification) the author set out to do a similar analysis 
for P5. It is by means of this type of analysis that R5 was quickly obtained. 

The errors found from chapter nine were found significantly later. Because 
the errors made with these equations are consistent with the previous errors, 
little analysis has been done other than identifying the errors. 

The author finds it convincing that Poincare could have made an error of 
omission on his fourth example. The relation between the third and fourth 
examples strongly supports this belief and the fact that of many plausible mod- 
ifications of Poincare's fourth example, R4 is the simplest makes it convincing. 
However, the fifth example does not fit within this genre of system and it is 
not clear that the error made was one of omission. The further errors in chap- 
ter nine make it convincing that the remaining errors where merely standard 
mathematical errors. 

The author would like to conclude with a quotation of Poincare: 

"How is an error possible in mathematics? A sane mind should not be guilty 
of a logical fallacy, yet there are some very fine minds incapable of following 
mathematical demonstrations. Need we add that mathematicians themselves 
are not infallible?" pp 
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